Double-stage continuous-discontinuous superconducting phase transition in the Pauh 
paramagnetic limit of a 3D superconductor: the URu2Si2 case 
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The sharp suppression of the de-Haas van-Alphen oscillations observed in the mixed superconduct- 
ing (SC) state of the heavy fermion compound URu2Si2 is shown to confirm a theoretical prediction 
of a narrow double-stage SC phase transition, smeared by fiuctuations, in a 3D paramagnetically- 
limitted superconductor. The predicted scenario of a second order transition to a nonuniform 
(FFLO) state followed by a first order transition to a uniform SC state, obtained by using a non- 
perturbative approach, is also found to be consistent with recent thermal conductivity measurements 
performed on this material. 



The competition between orbital and spin pair- 
breaking in strongly type-II superconductors in the Pauli 
paramagnetic limit is known to control the occurrence of 
discontinuous SC transitions 0,13] at sufficiently low tem- 
peratures and high magnetic fields. It was found recently, 
using perturbation expansion in the SC order parameter 
0, that in a clean 3D system, the normal-to-SC phase 
transitions at low temperatures are of second order, with 
a SC phase spatially modulated along the field direction 
0, 0], whereas the transition line from nonuniform-to- 
uniform SC state was found to be of the first order. This 
conclusion was reached, however, on the basis of pertur- 
bation theory, which might not be valid under the present 
circumstances due to the following reasons: (1) The jump 
of the SC order parameter to a finite value at the first or- 
der phase transition, and (2) the oscillatory dependence 
of the quartic and higher order terms in the expansion 
on the modulation wave number, which makes the uti- 
lization of a uniquely defined expression for the SC free 
energy meaningless within perturbation theory. 

The heavy fermion superconductor URu2Si2, whose 
Fermi surface (FS) may be characterized as 3D [6|, pos- 
sesses characteristic FS parameters which favor strong 
spin pair breaking. In this material a sharp rise of the 
thermal conductivity with the decreasing magnetic field 
just below Hc2 at low temperatures was reported very 
recently f?], indicating the existence of a jump in its elec- 
tronic entropy associated with a first-order phase transi- 
tion. Furthermore, earlier magneto-oscillations measure- 
ments on this material ^ revealed a very sharp damping 
of the de Haas-van Alpen (dHvA) effect just below Hc2 
which seems to correlate with the anomaly observed in 
the thermal conductivity. 

In this communication we present results of a non- 
perturbative approach, which establishes the sharp, 
double-stage transition picture, conjectures in Ref.[3|], 
and argue by means of a detailed theoretical analysis of 
the experimental dHvA data, that the predicted double- 
stage transition is realized in URu2Si2. The proposed 
model is also shown to be consistent with the anomaly 
in the thermal conductivity data reported in Ref.[7|. 

We start by writing an expansion of the thermody- 



namical potential (TP), 51, in the SC order parameter, 
A (r), using BCS theory for an isotropic 3D electron gas 
with the usual s-wave electron pairing, as presented in 
3]. The use of conventional pairing was made for the 
sake of simplicity. This is justified in the clean limit 
considered here since the relevant results have shown in 
Ref. ^3] to be independent of the type of electron pairing. 
Thus we write: 



Qint 



^2n (Ao) : 



^2n ^ / d3{r}r2„({r},Ao)if2„({r}) 



(1) 



where gint is the effective BCS coupling constant, V is 
the volume: 

r2„({r} , Ao) = ff*(ri,r2)5(r2,r3)...5r*(r2„_i,r2„) 
xg(r2„,ri)A(ri)A*(r2)...A(r2„_i)A*(r2„) (2) 



and: 



K2n{{A) = fcsT'2^Goj^(ri,r2,t^,.)GoT(r2,r3,cj^) 

■•■Goi (r2Ti-i, r2«, t^z.)Got (r2„, ri, w^) (3) 

Here A^ ^ V'^ J d^r, |A(r,)|', and {r} -{n, r2„} 
denotes the entire set of position vectors for a cluster con- 
sisting of n electron pairs. Note that for convenience we 
incorporated the gauge factors, g(ri, r^+i), of the Green's 
functions, Go|x (r^, r^-i-i, Wi^), for a free electron in a uni- 
form magnetic field, into the vertex part, r2„, so that the 
effective kernel i^2n is given in Eq. ([3]) by a product of the 
gauge invariant Green's functions, Go|j,(ri, r^+i, cj^). A 
useful expression for such a Green's function for a posi- 
tive Matsubara frequency, ly > 0, can be written as: 



GoTi(ri,r2,w^.) 



2-K 



dre 



exp 



2(1 



2 



where pi,2 = r_L2 — rj_i , with rj_i, r_L2-the projections 
of the initial and final electron position vectors, respec- 
tively, on the (x—y) plane perpendicular to the magnetic 
field. This expression is obtained after summation over 
the Landau level (LL) index n — 0,1,..., of the single- 
particle energy, Sn^^Jhwc — n + k'^ ^ g — np + itu^ for a 
spin up ( or down ) electron in a magnetic field H = i/z 
, with a cyclotron frequency Uc = eH/m*c , Zeeman spin 
energy ^eH/moc, and g- factor g = m* /mo, with m*, and 
toq- the effective mass and free electron mass respectively. 
Here /c^ = h?k'^/2m*hujc, np = jjb/hwc, /x-the chemical 
potential (« Ep -Fermi energy), and = cj^/wcwith 
= nksT {2i^ + 1) /h. For negative Matsubara fre- 
quencies, < 0, a similar expression can be derived by re- 
placing T with — T. In what follows we will express space 
coordinates and momenta in units of an — \J ch/eH 



and 



''H 



respectively. The SC order parameter is as- 



sumed to take the form, A(r, z) = Amaxe"'^V'o(a^, y), 

/ \ 1/2 

where A^-.^^^ = (a^j ' ^'"^"^ fo{x,y) describes (in 

the symmetric gauge) an hexagonal vortex lattice with 

inter- vortex distance, Qx = -0^- Here e"'^ is a Fulde- 
Ferrel (FF) Q modulation function along the magnetic 
field direction, controlled by the wave-number, q. 

The vertex part r2„({r},Ao), Eq.(l2]), is a violently 
oscillating function of the lateral relative electronic co- 
ordinates, which interferes strongly with the oscillatory 
electronic kernel i<'2n({r}), Eq.([3]). Multiple integration 
over these coordinates yields gross cancellations except 
near stationary configurations, which restrict all 2n elec- 
tronic position vectors to a relative proximity region of 
size of a magnetic length Other contributions 

to this integral, arising from non-stationary, separately 
paired configurations, become increasingly important in 
random vortex lattices where the phase coherence re- 
sponsible for the constructively interfering configurations 
breakdown 0. Whereas the small oscillatory (high har- 
monic in 1/7?) part of the TP is strongly infiuenced by 
these non-local contributions, their influence on the much 
larger non-oscillatory (zero harmonic in \/H ) compo- 
nent is not important (see Ref.Q). The great advan- 
tage of using the local approximation in Eq.([T]) is in 
its factorization with respect to the relative coordinates 
and its apparent independence of the center of mass co- 
ordinates, which enable us rewriting the integrand in 
Eq. ([T]) as a separable product of effective single electron 
Green's functions. The corresponding n-th order term, 
^2n , can be thus written as a 3D integral over the center 
of mass momentum in an effective two-particle Green's 
function, raised to the n-th power, by performing an ap- 
propriate Fourier-transformation, namely: fl2n (Amax) = 

fcgT Qx a2w t a2 
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where A^^^^ 
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and: 
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^^pi[a>.(k,fc,|5,g)$:(k,fc. 



Jo 



)|k|= 



with: Ck.fc, (5,9) = ^ |k| + 5 (^^ ~ 2 ) - 9 - np , 
k — xkx + yky . The resulting perturbation series can 
be easily summed to all orders, provided the reduction 
pre-factor , arising from the overlap integral of n 
lowest LL orbitals is represented as a Gaussian inte- 
gral: Jp°° exp (— nu^) du . In the quasi-classical limit 

T <C 1, the Gaussian approximation e^2(i+''^~'= )|k| ^ 

1 2 1. |2 

e~i'^ I'^i , accounts for the diamagnetic pair-breaking 
whereas all quantum corrections (including quantum 
magnetic oscillations), which arise near the lattice points 
T = 2ttI ,1^0, are neglected. 

It is convenient to normalize all energies by TrksTco 
, where Tco is the transition temperature at if = , 



so that: A 

1^,2 _ n^\k\^/2ra' 



77 _ £L 



max — TTkeT^o ' ^ ~ irfeeTeo ' 
■2 _ h^kl/2m 



h-'q^ /2m* 



, k^ — 



and g 



TTkBTaO ' 



ujc20 = £Hc2q/iti*c, and i?c20 is the theoretical upper 
critical field at T = (in the absence of spin splitting). 
Performing the integration over r, the resulting expres- 
sion for TP can be written in the form: 
-2 



n 



A 



TTfl 



kdkdkz In ( 1 



(■ 



e"" A 

^ max 1/ z> 



(4) 



(2cJc2oo) k 



t{2v +l)±i[k^ + {kz±qf ±gb-Ji 



{2ujc2ob) 

where b = H/Hc20, t = T/Tco , u Jc2q = hhJc2o/'^kBTco , 
(/ i^))u ^ r duf («) , i = V^E.>o 2^1, and 



So = 0.137 [kp^ {0)f with C (0) = 0.18hvp/kBT, 



cO 



the zero temperature coherence length. 

In the limit of zero spin splitting and in the absence of 

FF modulation {g ^ q ~ 0) the effective pairing parame- 
2 

ter, Xi, (fc, kz) = Aj^j^^$+$j7, is always real and positive. 
Under these circumstances the general form of ft ^A^^^^) 
at 77 < Hc2o, as expressed in Eq.(l4]), possesses the single 
minimum structure characterizing the usual GL theory. 
For g ^ , {k,kz) can become complex (a feature that 
can be "healed" by the presence of the FF modulation 



^2 



wavenumber g ), so that the general form of Q 

2 

may show a maximum at small Aj^g^^ which is followed 
2 

by a minimum at larger A^^^. The initial maximum re- 
flects the competition between the increasing spin param- 
agnetic energy and decreasing SC pair-correlation energy 
as the number of spin-singlet Cooper-pairs is increased. 

Typical results of the calculated TP using Eq.Q for 
non-zero spin-splitting are show in FigjT] As discussed 



3 




FIG. 1: (color online) n vs. Amax (in units of S.qV) for a 
uniform ( g = ) SC order parameter (solid curves) and 
for the corresponding FF modulated order parameter (dashed 
curves) at various magnetic field values b near the SC tran- 
sition at temperature t — 0.02. The selected g- factor is 
g = 1.8. Note the second-order transition to the FF state 
at 6 « 0.195 and the first-order transition to the uniform SC 
state at b ~ 0.185. 
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FIG. 2: (color online) Field dependence of the self-consistent 
order-parameter amplitude Amax at a temperature t = 0.02 , 
well below the tri-critical point {ttric ~ 0.4), for g = 1.8. 



above, the restriction to SC states with q — , repre- 
sented in Figdlby the solid curves, leads to formation of 
a maximum at small Amax and a local minimum at larger 
Amax by the strong spin splitting effect as the magnetic 
field is reduced. Upon further field decrease the minimum 
becomes global and so should drive a first order normal- 
to-SC phase transition. Allowing for states with q 0, 
however, the unusual feature (i.e. Im X^, (k, kz) 7^ ) as- 
sociated with the strong spin splitting effect is " healed" , 



and the usual single minimum picture is restored (see the 
dashed curves in FiglJ). Thus, instead of the "expected" 
first-order transition to a uniform SC state one finds a 
second-order phase transition to a nonuniform (FF) SC 
state. 

However, due to its compensation effect, the FF mod- 
ulation significantly reduces the equilibrium SC free en- 
ergy with respect to its uniform counterpart (compare 
the dashed curves to the corresponding solid ones). As 
a result, the field range of stability of the modulated 
phase is quite small. Thus, by slightly reducing the 
field below the second order normal-to-SC transition the 
(7 = state becomes energetically more favorable and 
the system transforms from the nonuniform to a uni- 
form SC state via a first-order phase transition. One 

should note that the second order perturbation theory 

2 1 4 

with fi/FSg = a {q) A^ax + ^max: quantita- 

tively correct only for small values of the order param- 
eter, namely for Amax ^ -Ij and, therefore, cannot be 
applied to the first order transition. 

2 

Fig[2] shows the calculated self-consistent Amax (which 
is optimized with respect to g ) as a function of magnetic 
field at a temperature well below the tri-critical temper- 
ature and for characteristic parameters corresponding to 
URu2Si2. The initial build-up of the SC order parame- 
ter in a narrow region following the second-order phase 
transition and the pronounced jump at the next (first 
order) transition are apparent. The overall width of the 
two-stage transition {Ab/b « 0.06 ) is in good agreement 
with the width of the sharp structure observed experi- 
mentally in the thermal transport measurements 

FigO exhibits the result of a detailed fitting proce- 
dure of our calculation to the experimental data of the 
dHvA oscillations observed by Ohkuni et. al [(| in 
the mixed state of URu2Si2. The data shows a very 
sharp reduction in the amplitude of the dHvA oscilla- 
tion, Asc I just below H = 2.78 T , similar to the 
step-like structure observed in the thermal conductiv- 
ity measurements 0- For fitting the measured relative 
signal, Asc/ An ( being the theoretical dHvA ampli- 
tude as extrapolated from the normal state) in the region 
above the sharp damping interval, we exploit the fluctu- 
ating vortex lattice model described in Refs.0,[l3] and 
write: ln(AscM„) = -^-^g^ (A^) , where (A^) 
is the mean-square order parameter for the modulated 
FF state in the vicinity of the second-order phase tran- 
sition,i.e.: (A^) = -^A^^^ (l + ^1 + u^/x^) , with 

X = -j= \/ Amax; and = 0.51 for a 3D system 
(see Ref.[io|). Here (3 is the coefficient of the quartic 
term obtained in the expansion of in Amax near 
the SC transition, and Vq = n'^ajj^^/kp- The selected 
value of Ep has been determined from the experimen- 
tally observed, dominant dHvA frequency, Fa ~ 10^ T, 
corresponding to the nearly spherical band 17- hole Fermi 
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FIG. 3: (color online) Logarithm of the dHvA amplitude ra- 
tio Asc/^n reported in Ref.^Q] (full circles) and the values 
of — y(A^) , y = as obtained from our calcula- 

tion near the second order phase transition (solid line), as 
functions of 1/H. The dashed broken straight line is our 
mean-field result for — ^A^^^ around the first-order transi- 
tion. The doted line represents the extrapolation of — x (A^) 
for the modulated FF state to the low-field regime. 



surface reported in Ref.[6|. 

2 

The values of Aj^g^^ at H < Hc2 were determined in 
this calculation by minimizing the TP with respect 

2 

to both Ajjjjjx and q (see Figl2]), whereas aX H > Hc2 
they were obtained by the analytical continuation of 

2 

the perturbative mean field expression, A^^^^^ = —a/ (3 
into the region —a/ (3 < 0, where a and /3 are the 
quadratic and quartic coefficients respectively in the per- 
turbation expansion of il ^^A^^^^^ . The optimal value 
of q in this region was selected by minimizing with re- 
spect to q the free energy obtained from the functional 
integral of e~^(^"'»'')^'^^"^ over A^^^^, which amounts in 
the Gaussian approximation to minimizing a with re- 
spect to q. In the fitting procedure we have exploited 
the interpolation formula (A^) ~ ^D"^ (l — H / Hc2) + 

^Joi] + jD^ (1 - H/Hc2f^ with the best fitting ad- 
justable parameters Dq — 4.2 K, Hc2 — 2.76 T , and 
ai = .06 K (see Figl3]). The resulting value of the order 
parameter is in a reasonably good agreement with the 
zero-field gap parameter Abcs = 1.76kBTcQ — 2.5 K 
obtained within BCS theory with the experimental value 
of r,o (= lAK). 

Using the resulting parameters and the field depen- 

2 

dent Ajjjjj^ shown in FiglSj we have calculated the jump of 
In (Asc/^n) at the first order phase transition within our 
nonperturbative mean field theory. The agreement with 
the experimental data is good, keeping in mind that ther- 
mal fluctuations and inhomogeneous broadening could be 
responsible for the observed smearing of the theoretical 



discontinuous transition. Nevertheless, the sharp down- 
ward deviation of the experimental data in Figl3] from 
the extrapolation of the calculated — x (A^) for the mod- 
ulated FF state to the low-field regime provides a strong 
evidence for the two-stage nature of the SC transition. It 
should be stressed that the relative size of the jump in 
Amax obtained in our calculation is nearly independent 
of the various parameters involved, provided the temper- 
ature t is well below the tri-critical point ttric- 

In conclusion, using a non-perturbative approach, we 
have firmly established our early conjecture concerning 
the SC transition in a 3D strongly type-II superconduc- 
tor in the paramagnetic limit, and show that the dHvA 
effect observed in the mixed SC state of URu2Si2f^ pro- 
vides a clear experimental evidence for the double-stage 
nature of this transition, which is smeared by significant 
SC fluctuations effect. This flnding is consistent with 
the interpretation of a first-order phase transition given 
in Ref. 0| to the step-like structure observed in the ther- 
mal transport data of this material. We note that the 
unusual sign of the observed jump in the thermal conduc- 
tivity, which could be due to some peculiar quasi-particle 
scattering mechanism llj, is irrelevant to our main ar- 



gument, which associates this jump, irrespective of its 
direction, to the jump of the SC order parameter at the 
predicted first-order transition. 
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